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, Recent spectroscopic observations of planetary nebulae (PNe) in several ellipti- 

cal galaxies have revealed structural and kinematical properties of the outer stellar 
halo regions. In order to elucidate the origin of the properties of these planetary neb- 
ula systems (PNSs), we consider the merger scenario in which an elliptical galaxy is 
formed by merging of spiral galaxies. Using numerical simulations, we particularly in- 
1^ \ vestigate radial profiles of projected PNe number densities, rotational velocities, and 

(^jQ ■ velocity dispersions of PNSs extending to the outer halo regions of elliptical galax- 

' ies formed from major and unequal-mass merging. We find that the radial profiles of 

\^ , the project number densities can be fitted to the power-law and the mean number 

' density in the outer halos of the ellipticals can be more than an order of magnitude 

I higher than that of the original spiral's halo. The PNSs are found to show a signif- 

, icant amount of rotation [V/a > 0.5) in the outer halo regions {R > 5i?o) of the 

I ' ellipticals. Two-dimensional velocity fields of PNSs are derived from the simulations 

O I and their dependences on model parameters of galaxy merging are discussed in de- 

^ . tail. We compare the simulated kinematics of PNSs with that of the PNS observed in 

^ ' NGC 5128 and thereby discuss advantages and disadvantages of the merger model in 

. . , explaining the observed kinematics of the PNS. We also find that the kinematics of 

■ PNSs in elliptical galaxies are quite diverse depending on the orbital configurations 
k> ' of galaxy merging, the mass ratio of merger progenitor spirals, and the viewing angle 

. of the galaxies. This variation translates directly into possible biases by a factor of 

' two in observational mass estimation. However, the biases in the total mass estimates 

can be even larger. The best case systems viewed edge-on can appear to have masses 
lower than their true mass by a factor of 5, which suggests that current observational 
studies on PNe kinematics of elliptical galaxies can significantly underestimate their 
real masses. 
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ies individual (NGC 5128) — galaxies:kinematics and dynamics — galaxies:stellar 
content 



1 INTRODUCTION 

Structural, kinematical, and chemical properties of stars 
in the outer regions of galaxies, their halos, provide vi- 
tal clues to their global formation and evolutionary histo- 
ries. Although the physical properties of stellar halos have 
been extensively investigated mostly for spiral (e.g., the 
Galaxy and the M31) and dwarf galaxies in the Local Group 
(e.g., Eggen, Lynden-Bell, & Sandage 1962; Mould & Kris- 
tian 1986; Norris 1986; Freeman 1987; Durrell, Harris, & 



Pritchet 1994, 2001; Pritchet & van den Bergh 1998; Re- 
itzel, Guhathakurta, & Gould 1998; Grillmair et al. 1996; 
Chiba & Yoshii 1998; Davidge 2002), recent observations 
have started to reveal physical properties of stellar halos in 
giant elliptical galaxies, all of which are beyond the Local 
Group (e.g., Soria et al. 1996; Harris, Harris, & Poole 1999 
hereafter referred to as HHP; Harris & Harris 2000, HHQQ; 
Harris & Harris 2002, HHQ2; Marleau et al. 2000; Rejkuba 
et al. 2002; Gregg et al. 2004). These data prompt us to in- 
vestigate how the properties of stellar halos can be used to 
understand the formation of giant elliptical galaxies. 
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Beyond projected radii of ^ 2Rc, it becomes very difh- 
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cult to study the integrated properties of stars in ellipticals 
because of the low stellar surface brightness. Hence, stud- 
ies of stellar halos tend to rely on resolving individual stars 
or star clusters. While there are significant observational 
difficulties in revealing the three dimensional structure and 
kinematics of halo red giants and AGB stars, it has been pos- 
sible to use color-magnitude diagrams to derive their metal- 
licity distribution function (hereafter MDF) and constrain 
the past star formation and chemical evolution histories of 
their host elliptical galaxies (e.g., HHP; HHOO, HH02). Kine- 
matic properties of these stars, however, which can contain 
vital clues to both structure and formation of ellipticals, are 
still extremely difficult to obtain. Furthermore, these pre- 
vious observations have difficulties in deriving global stellar 
halo distributions, which are also considered to have valu- 
able information on the past merging histories of galaxies 
(e.g., fossil tidal streams), mainly because the halo fields 
investigated cover only a minor portion of the entire halo. 

Planetary nebulae (PNo) are a valuable complement to 
integrated light and red giant studies. They are powerful 
tools for studying the dynamical states of elliptical galax- 
ies, because we can more readily identify PNe through their 
bright [O HI] omissions linos and typically measure their ra- 
dial velocities to an accuracy of 15 km s~^ (e.g., Peng et 
al. 2004a, PFF04a) . Spectroscopic observations of the radial 
velocities of PNe in several neaxby elliptical galaxies there- 
fore have succeeded in deriving global mass distributions and 
velocity fields of those galaxies (CiarduUo et al. 1993; Hui 
et al. 1995; Arnaboldi et al. 1998; Mendez et al. 2001; Ro- 
manowsky et al. 2003; Napolitano et al. 2004; PFF04a). For 
example, kinematic information of PNe in the stellar halo 
of NGC 5128 at radii up to ~ 80 kpc allows the authors to 
derive the zero velocity curve (ZVC) in the two-dimensional 
velocity field of the PNS and thereby to discuss the triaxial 
shape of the mass distribution of this galaxy (PFF04a). This 
observed diversity in kinematics of PNSs is providing fresh 
clues to the origin of early-type galaxies. (Romanowsky et 
al. 2003; PFF04a). 

In spite of this importance of PNe studies in ellipti- 
cal galaxies, only a few theoretical and numerical attempts 
have boon made to discuss dynamical properties of PNSs in 
elliptical galEixies. For example, using kinematical data up 
to ~ 6Re for an heterogeneous sample of elliptical galax- 
ies, Napolitano et al. (2004) derived the dependences of the 
radial gradients of the mass-to-light-ratios on the B— band 
magnitude of the galaxies and thereby discussed whether 
the derived dependences can be consistent with elliptical 
galaxy formation models based on a ACDM model. How- 
ever, PNSs provide a very rich data set from which to study 
ellipticals. In particular, we are interested in exploring the 
2D distributions of structural and kinematical properties of 
PNSs, which can provide some vital clues both to the triax- 
ial shapes of the global mass distributions of galaxies and to 
their formation processes. The origin of 2D structural and 
kinematical properties of PNSs extending to the outer halo 
regions of elliptical galaxies remains unclear. Furthermore 
providing theoretical predictions on dynamical properties of 
PNSs will help to interpret the properties of PNSs that will 
be obtained in future systematic observations for elliptical 
galaxies. 

The purpose of this paper is thus to investigate struc- 
tural and kinematical properties of PNSs of elliptical galax- 



ies based on numerical simulations. In order to elucidate the 
origin of the observed properties of the PNSs, we adopt the 
merger scenario (Toomre 1977) in which elliptical galaxies 
are proposed to be formed by major merging of two spiral 
galaxies. The present numerical investigation is two-fold as 
follows. We first describes radial profiles of structural and 
kinematical properties of PNSs and two-dimensional (2D) 
velocity fields of PNSs in elliptical galaxies and their depen- 
dences on model parameters. We then compare the simu- 
lated kinematics of the PNSs with the corresponding obser- 
vations for NGC 5128 and thereby try to provide the best 
model for the PNS in NGC 5128 and discuss advantages 
and disadvantages of the model in explaining the kinemat- 
ics of the PNS sclf-consistcntly. The essential reason for our 
choice of the NGC 5128 PNS is that PFF04a have recently 
investigated radial velocities of the NGC 5128's 780 PNe 
(among 1141 PNe), which represents the largest kinemat- 
ical study of an elliptical galaxy to date and thus can be 
compared with our simulations in the most self-consistent 
manner. Furthermore, the present study is complimentary 
to those by Bekki et al. (2003, BHH03) and Beasley et al. 
(2003) which numerically and semi-analytically investigated 
the MDF of the stellar halo of NGC 5128 but did not inves- 
tigate the dynamical properties. 

In the present paper, we mainly investigate structural 
and kinematic properties of PNSs for a large radial extent in 
elliptical galaxies (0 ^ i? ^ 10-Rc), which include the outer 
faint stellar halos and the main bodies of the galaxies. Al- 
though previous numerical simulations investigated spatial 
distributions and the line-of-sight velocity distributions of 
stars in elliptical and SO galaxies formed from merging for 
R < 2.5Rc (e.g., Bendo & Barnes 2000; Cretten et al. 2001; 
Naab & Trujillo 2005), they did not investigate the dynami- 
cal properties for the outer halo regions of elliptical galaxies. 
Therefore, the present numerical results on the dynamical 
properties of the outer stellar halos including PNe {5Rc ^ R 
^ lO-Ro) may provide new clues to elliptical galaxy forma- 
tion. Physical properties of these halo stars may well have 
fossil information on angular momentum redistribution of 
stars (or conversion from galactic orbital angular momen- 
tum into intrinsic angular momentum) which is an essential 
physical process of galaxy merging. We therefore can expect 
that dynamical properties of stellar halos, which can be pro- 
vided by studies of PNSs, may depend strongly on physics of 
galaxy merging. Through the course of this paper, we thus 
interpret our simulations in the context of the capabilities of 
present-day PN surveys in order to facilitate the comparison 
between theory and observations. 

The plan of the paper is as follows: In the next sec- 
tion, we describe our numerical models for the formation of 
PNSs extending to the outer halo regions in galaxy merg- 
ers. In §3, we present the numerical results mainly on the 
final 2D distributions of structural and kinematical proper- 
ties in merger remnants (i.e., elliptical galaxies) for variously 
different merger models. In §4, we compare the simulated re- 
sults of PNSs with observations of the PNS in NGC 5128 
(Gen A). In §5, we discuss whether physical properties of 
stellar halos can give any constraints on galaxy formation, 
based on the present numerical results. In this section, wc 
also discuss possible different properties of PNSs in galaxies 
with different Hubble types. We summarize our conclusions 
in §6. 
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Table 1. Model parameters 



Model no. 


Orbital type 


Ml 




M2 


HI 


M3 


PP 


M4 


RR 


M5 


LA 


M6 


BO 


M7 


HI 


M8 


HI 


M9 


PO 


MIO 





2 MODEL 

2.1 Merger models 

We investigate the dynamical evolution of fully self- 
gravitating galaxies composed of stars and dark matter via 
coUisionless numerical simulations carried out on a GRAPE 
board (Sugimoto et al. 1990). Since our numerical methods 
for modeling dynamical evolution of galaxy mergers have 
already been described by Bekki & Shioya (1998) and by 
BHH03, we give only a brief review here. The total mass 
and the size of a disk of the merger progenitor spiral are 
Md and Rd, respectively. Henceforth, all masses and lengths 
are measured in units of A4d and Rd, respectively, unless 
specified. Velocity and time are measured in units oi v = 
{GMd/Rdy^ and tdyn = [Rd/GMdY^^, respectively, where 
G is the gravitational constant and assumed to be 1.0 in the 
present study. If we adopt Md = 6.0 x 10^° M© and _Rd = 
17.5 kpc as a fiducial value, then v — 1.21 x 10"^ km/s and 
^dyn = 1.41 X 10* yr, respectively. The merger progenitor 
spiral is composed of a dark matter halo, a stellar disk, a 
stellar bulge, and a stellar halo. 

The mass ratio of the dark matter halo to the stellar 
disk in a disk model is fixed at 10 for all models. We adopt 
the density distribution of the NFW halo (Navaxro, Frenk 
& White 1996) suggested from CDM simulations: 

where r, po, and rs are the spherical radius, the central den- 
sity of a dark halo, and the scale length of the halo, respec- 
tively. The value of rs (0.6 in our imits for c = 10) is chosen 
such that the rotation curve of a disk is reasonably con- 
sistent with observations. The mass fraction and the scale 
length of the stellar bulge are fixed at 0.17 (i.e., 17 % of the 
stellar disk) and 0.04 (i.e., 20 %) of the scale length of the 
stellar disk), respectively, which are consistent with those 
of the bulge model of the Galaxy. The initial total mass of 
the stellar halo is 0.01 in our units for all models and has a 
number density distribution (p(r) ~ r~^'^) like that of the 
stellar halo of the Milky Way (e.g., Chiba & Beers 2000). 

The total number of coUisionless particles in each sim- 
ulation is 112680 for pair mergers and 281700 for multiple 
merger (described later). The total particle number for dark 
matter halo, stellar disk, bulge, stellar halo, and old halo 
GCs are 30000, 20000, 3340, 2000, and 1000, respectively, 
for a disk. Gravitational softening lengths in the GRAPE 
simulations are set to be fixed at 0.15 (in our units) for the 
dark matter, 0.02 for the stellar disk and bulge, 0.09 for the 
stellar halo for each simulation. Since we mainly analyze 



m2 Comments 
- isolated disk 

1.0 fiducial model 
1.0 
1.0 
1.0 
1.0 
0.1 
0.3 

0.5 NGC 5128 model 
1.0 Multiple merger 

structural and kinematical properties of PNSs for the scale 
of 1 — lOOkpc in galaxies, the resolution of the simulations 
(~ 350 pc) is quite reasonable. 

We investigate both pair mergers with the number of 
the progenitor spirals (A'g) equal to 2 and multiple ones 
with A'g = 5. Since our main interest here is the formation 
of elliptical galaxies by pair mergers, the details of the mod- 
els and the results for the multiple mergers are given in the 
Appendix B. In the simulation of a pair merger model, the 
orbit of the two spirals is set to be initially in the x-y plane 
and the distance between the center of mass of the two spi- 
rals (rin) is 6iid (105 kpc). The pericentre distance (rp) and 
the orbital eccentricity (ep) arc assumed to be free param- 
eters which control orbital angular momentum and energy 
of the merging galaxies. For most merger models, rp and Cp 
are set to be 1.0 (in our units) and 1.0, respectively. The 
spin of each galaxy in a merging and galaxy is specified by 
two angles 9i and (f>i (in units of degrees) , where the suffix i 
is used to identify each galaxy. Here, 9i is the angle between 
the z-axis and the vector of the angular momentum of the 
disk, and is the azimuthal angle measured from x axis to 
the projection of the angular momentum vector of the disk 
onto the x-y plane. 

We specifically show the following five different and rep- 
resentative models with different disk inclinations with re- 
spect to the orbital plane: A prograde-prograde model rep- 
resented by "PP" with Oi ^0,02= 30, </-i = 0, and (?!-2 = 0, a 
retrograde-retrograde ("RR") with 6i = 180, d-z = 150, 01 = 
0, and (f>2 = 0, and a highly inclined model ( "HI" ) with di = 
30, 6>2 = 120, 01 = 90, and 02 = 180. The model with lower 
orbital angular momentum (rp = 0.2) and a HI orbital con- 
figuration is labeled as "LA" whereas the model with bound 
orbit with ep of 0.7 and a HI orbital configuration is labeled 
as "BO". 

The total masses of elliptical galaxies can be underes- 
timated, if kinematical properties of PNSs are used for the 
mass estimation (e.g., Dekel et al. 2005). We also find that 
the masses of ellipticals can be significantly underestimated 
if radial velocity profiles of PNSs are used under the assump- 
tion of the virial theorem. We just briefly discuss this point 
in the discussion section §5.3, because our main focus is on 
2D density and velocity fields of PNSs in elliptical galaxies. 



2.2 Main points of analysis 

We mainly investigate 2D density and velocity fields of PNSs 
and compare the results with the corresponding observations 
in a fully self-consistent manner. We accordingly use Gaus- 
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Figure 1. The two-dimensional (2D) PNe density field of the 
stellar halo in the isolated disk model Ml projected onto the X-Y 
plane. The stellar halo region with the size of 17.5 kpe is divided 
into 40 X 40 eells and the projected PNe density at each cell is 
estimated by assuming cipNe of 9.4 X 10~^ PNe per Lq in B-band 
(CiarduUo et al. 1989). The original spherical distribution of the 
PNS is clearly seen in this smoothed 2D density field. 



sian smoothing similar to that adopted in observations (e.g., 
PFF04a) in order to derive smoothed density and velocity 
fields of the simulated PNSs. The details of the methods to 
derive the smoothed fields are given in the Appendix A. Fig- 
ure 1 shows how the smoothed PNe density field looks like 
for the stellar halo consisting only of 1000 discrete stars in 
an isolated disk model Ml (viewed from face-on). 

Although the MDFs of stellar halos (thus PNe) can give 
some constraints on any theory of elliptical galaxies forma- 
tion (e.g., HHOO; HH02; BHH03), we do not intend to discuss 
them in this paper. This is firstly because (1) the importance 
of the MDFs have been already discussed in previous papers 
(e.g., BHH03) and secondly because there is no currently 
available data set of the MDFs of PNe for nearby giant el- 
liptical galaxies. Instead, we mainly investigate (1) radial 
profiles of projected PNe number densities, (2) 2D fields of 
line-of-sight velocity (Vios) and velocity dispersion (cr), (3) 
rotation curve profiles (Vrot), and (4) radial profiles of a and 
Kot/o", for the PNSs of merger remnants. Although we in- 
vestigated the above four points for 32 merger models with 
different m2 and orbital configurations, we show the results 
of 9 representative models among them. Table 2 summarizes 
the model parameters used for these 9 models. We show the 
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Figure 2. The mass distributions of stars projected onto the 



X-Y plane (face-on view, left) and the X-Z one (edge-on view, 
right) for the fiducial model M2 at T = 4.5 Gyr. 



results of the models at T = 4.5 Gyr (at final time step), 
where the time T represents time that has elapsed since the 
merger progenitor disks begin to merge. 

We firstly show the results of the "fiducial model", 
model M2, which show typical behaviors of stellar halo (thus 
PNS) formation in major galaxy merging (§3.1). The results 
of this model can be regarded as generic ones for elliptical 
galaxies formed from major merging. Secondly we show the 
parameter dependences of the results in §3.2 with special 
emphasis on the dependence of projected PNe number den- 
sities and the 2D velocity fields on orbital configuration and 
7712. Thirdly, we present the best model which can repro- 
duce the most reasonably well a number of dynamical prop- 
erties of the PNS observed in NGC 5128 (§4.1). It should be 
stressed here that the best model is the best model among 
the present 32 models, thus is not the one which reproduce all 
of the observed properties of the PNS fully self-consistently. 
We try to understand advantages and disadvantages of the 
present best model by comparing it with observations, and 
discuss whether (and how) physical effects that are not in- 
cluded in the present models (e.g., star formation) can be 
important for more successful reproduction of the observa- 
tions. The best model is model M9 with ep = 1.0, rp = 0.05 
in our units, 9i = 0, ^2 = 80, 01 = 0, and 02 = 0, and 7712 
= 0.5 (The orbital configuration is hereafter referred to as 
"PO"). 



3 RESULTS 

3.1 The fiducial model 

3.1.1 Density fields 

Figure 2 shows the projected stellar mass distributions of 
the fiducial model M2 at T = 4.5 Gyr in which the final 
merger remnant shows of 5.2 kpc, Lb = 3.7 x 10^° Lq, 
and Np-Me = 350 within ~ 5-Re. The merger remnant does 
not show any clear signs of stellar substructures (i.e. tidal 
tails and plums) in its halo region for the two projections 
owing to the apparently completed dynamical relaxation by 
this time (T = 4.5 Gyr). It is also possible that the much 
less remarkable substructures are due to the adopted smaller 
particle number (an order of ~ 10^). The outer stellar halo 
{R > 5i?o) in the edge-on view is more spherical (e ~ 0.3) 
than the main body of the elliptical (e ~ 0.2 at i? ~ 2i?c). 



Density field (Face — on) Density field (Edge — on) 
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Figure 3. The 2D PNe density field for the face-on view (left) 
and the edge-on view (right) in the fiducial model at T = 4.5 Gyr. 
The 2D fields are produced based on the stellar mass distributions 
shown in Figure 2. 




R (kpc) 

Figure 4. Radial profiles of the projected PNe number densities 
(PpNc) for the fiducial merger model M2 (thick) and the isolated 
disk model Ml (thin). Note that the elliptical galaxy formed in 
M2 shows significantly higher ppNe in its outer region (R > 20 
kpc) compared with the spiral in Ml. 
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Figure 5. The 2D velocity fields for the face-on view (left) and 
for the edge-one one (right) for the fiducial model M2 at T = 4.5 
Gyr. 



Figure 6. The same as Figure 5 but for the 2D velocity disper- 
sion fields. 



The outer stellar halo is slightly more flattened than the 
dark matter halo with e = 0.15 at i? = 5Rc in the edge-on 
view and the major-axis of the stellar halo nearly coincides 
with that of the dark matter halo. The isophotal shape of 
this remnant formed from purely colhsionless major galaxy 
merging can show the negative sign of the 04 parameter 
(e.g., Bekki & Shioya 1997), and thus this remnant can be 
morphologically identified as a boxy elliptical galaxy. 

Figures 3 and 4 describe the 2D density fields of the 
PNS in the elliptical and the radial profile of the PNe dis- 
tribution, respectively. Owing to the adopted assumption 
of constant qpnc (i.e., constant luminosity-specific number 
densities of PNe), the 2D density fields of the PNS follows 
the stellar density ones of the elliptical and therefore shows 
a more fiattened shape in the edge-on view. The radial den- 
sity profile is best fitted to the power-law with the exponent 
of —2.7 for ^ R ^ 40 kpc (~ 8Re), which is slightly steeper 
than the observed slope of -2.5 for the PNS of NGC 5128 
(PFF04a). The projected PNe density of the elliptical in its 
outer stellar halo {R > 5-Rc) is more than two orders mag- 
nitudes higher than those of the merger progenitor spirals, 
though the PNe densities in the central regions of galaxies 
are not so different between the elliptical and the spirals. 
The formation of such a high density stellar halo is demon- 
strated to be closely associated with redistribution of disk 
stars through tidal stripping of the stars during major galaxy 
merging (BHH). 

Figure 5 shows the 2D fields of Vios (i-e., line-of-sight- 
velocity) of the merger remnant for ^ i? ^ 35 kpc (~ 7-Rc) 
in the face-on and edge-on views. Clearly, the PNS shows 
a significant amount of rotation (~ 100 km s~^) along its 
major axis (i.e., the X-axis) for the edge-on view. Given the 
fact that the stellar halos of the merger progenitor spirals 
are assumed to have no net rotation, the rotation in the 
outer halo of the elliptical can be considered to be obtained 
during angular momentum redistribution of major galaxy 
merging. This angular momentum transfer processes have 
been already well discussed in previous papers (e.g., Barnes 
1998 for a review). 

The axis of rotation in the edge-on view however does 
not coincide with the minor axis of the PNe density field 
shown in Figure 3: Radial velocity gradient can become max- 
imum if we measure it along a line connecting two points (X, 
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Figure 7. The radial profile of the rotation curve Vrot (upper) 
and the velocity dispersion <t (lower) for the dark matter halo 
(thick) and for stars (thin) in the fiducial model M2 at T = 4.5 
Gyr. Viot and a are estimated for the model projected onto the 
X-Y plane (i.e., the face-on view). The methods to derive Vrot 
and (T and their error bars are described in detail in the Appendix 
A. 

Z) — (—35, 35) and (35, —35) kpc for the edge-on 2D ve- 
locity field. Minor axis rotation can be clearly seen in the 
outer halo region of the elliptical for the edge-on 2D veloc- 
ity field and the velocity variation along the minor axis (i.e., 
the Z-axis) is complicated. Also it should be stressed that 
the edge-on 2D field has several regions with large absolute 
magnitudes of Vios (e.g., {X, Z) = (35, —8) and (—20, 3) 
kpc). Such inhomogeneity in the velocity space may indi- 
cate that dynamical relaxation is not completely ended in a 
real term even a few Gyr after the coalescence of two spirals 
in the fiducial model: Fossil evidences of past violent relax- 
ation can be proved in velocity structures of outer stellar 
halos in elliptical galaxies. 

The face-on 2D velocity field also shows rotation along 
the major and the minor axes of the projected density distri- 
bution of the PNS, though the rotation in the outer stellar 
halo for the facc-on 2D field is less remarkable compared 
with the edge-on one. The global appearance of the face-on 
2D velocity field is smoother than that of the edge-on one 
and complicated velocity variation along the minor axis of 
the PNS can not be clearly seen. Thus these results imply 
that (1) PNSs in elliptical galaxies can show a significant 
amount of rotation in their outer halo regions, (2) minor- 
axis rotation of PNSs is one of important characteristics of 
elliptical galaxies, and (3) the 2D velocity fields can be dif- 
ferent for different viewing angles. 

As shown in Figure 6, both the facc-on and the edge- 
on 2D velocity dispersion (ct) fields show (1) higher velocity 
dispersion in the inner regions of the elliptical, (2) a shal- 
lower gradient of velocity dispersion along the major axis 
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Figure 8. The same as Figure 7 but for the edge-on view (i.e., 
the X-Z plane). 



(i.e., the X-axis) compared with that along the minor one, 
and (3) globally flattened shapes, in particular, in outer stel- 
lar halo regions. These three can be regarded as important 
kinematical characteristics of stellar halos (PNSs) in ellipti- 
cal galaxies formed by major merging. The flattened velocity 
dispersion flelds are due to anisotropic velocity dispersion 
in the triaxial mass distribution of the elliptical galaxy in 
this fiducial model. The dispersion field appears to be more 
homogeneous in the face-on view than in the edge-on one, 
though some remarkable substructures that can result from 
yet incomplete dynamical relaxation of this system can be 
seen in both projections. The results shown in Figures 4, 
5, and 6 clearly suggest that even if the 2D density fields of 
PNSs appear to be quite regular and axisymmetric, their 2D 
velocity and velocity dispersion fields can be more inhomo- 
geneous and less axisymmetric: The kinematical properties 
of a stellar halo in an elliptical galaxy can provide some ev- 
idences of the past dynamical processes of its formation, if 
penetrative analysis of such "fossil record" properties can be 
done. 

Figures 7 and 8 show the radial profiles of rotational 
velocities (Kot) and velocity dispersions (ct) for the dark 
matter halo and the stellar components (i.e., PNe) in the 
elliptical of the fiducial model. It is clear from these figures 
that (1) the dark matter halo shows a larger velocity dis- 
persion (by a factor of 1.1 — 2.0, depending on radius) than 
the PNe, (2) the PNe shows a significant amount of rotation 
40 — 60 km s~^ in the entire halo region {R > 57?e) of the 
elliptical for the edge-on view, (3) the maximum Kot in the 
halo region is higher in the edge-on view than in the face-on 
one, (4) the radial gradient of Kot is quite different between 
different projections, and (5) the radial velocity dispersion 
profiles (e.g., the central a and the gradient) do not differ 
significantly with each other. The above result (1) strongly 
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Figure 9. The radial profiles of Vrot/c (upper) and ttie normal- 
ized rotation curve Vrot/o'o (lower) for the fiducial model viewed 
from 12 representative, and different viewing angles. Note that 
there is a trend of increasing K-ot/c and Vrot/fo with radius. 



suggests that wo can significantly underestimate the total 
mass of an elliptical galaxy, if we use only the velocity dis- 
persion data of the PNS and thereby derive the total mass 
based on the scalar virial theorem in which the mass is lin- 
early proportional to the product of -Re and a: We need to 
derive the total mass of an elliptical by using both dispersion 
and rotation data of the PNS. The result (2) is consistent 
with previous numerical simulations (e.g., Heyl et al. 1996). 

Given the fact that the radial dispersion profile of the 
PNS decreases with radius in the fiducial model, the above 
result (2) means that the PNS shows a higher value of Kot/o" 
in the outer halo of the elliptical, (e.g., Kot/o" of ~ 0.7 at 
R ~ 5Jie). This larger Kot/f also results from the redis- 
tribution of angular momentum of disks stars during major 
merging (i.e., conversion of initial orbital angular momen- 
tum of merger progenitor disks into final intrinsic one of 
the merger remnant), as explained for Figure 5. This result 
suggests that larger Kot /c in the outer stellar halo of an el- 
liptical galaxy can be a fossil evidence that the elliptical was 
formed from major galaxy merging. The derived PNS with 
larger Kot/o" in the outer halo region is in a striking contrast 
to the Galactic stellar halo with little rotation (e.g.. Free- 
man 1987), which implies that stellar halo kinematics can be 
significantly different between spirals and ellipticals. Recent 
numerical simulations based on the currently favored cold 
dark matter (CDM) theory of galaxy formation have shown 
that the Galactic stellar halo with little rotation is formed 
at high redshift {z > 1) by dissipative and dissipationless 
merging of smaller subgalactic clumps and their resultant 
tidal disruption in the course of gravitational contraction of 
the Galaxy (Bekki & Chiba 2000, 2001). Thus the above 
result (2) implies that differences in formation processes of 
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Figure 10. Vrot /<t aX R = 5Re as a function of T^ot /cr at R = 2Re 
(upper) and Vrot /fo at i? = bRe as a function of Kot /fo at R = 
2Rc (lower) for the fiducial model viewed from 12 representative, 

and different directions. 



stellar halos between different typos of galaxies (e.g., spirals 
vs ellipticals) can be rcfloctod on the stellar halo kinematics. 
The above results (3) - (5) imply that radial profiles of rota- 
tional velocities of PNSs in elliptical galaxies depend more 
strongly on viewing angles than those of velocity dispersion. 

Figure 9 describes how Vrot/cr and Kot/o"o depend on 
radius in the PNS of the fiducial model for 12 representative 
projections. There can be seen a trend of increasing Kot/c 
and Kot/cTo with increasing radius, though the radial pro- 
files rise and fall with radius significantly. The higher Kot/c 
and Vrot/co of the PNS in the outer halo region of the el- 
liptical for most of the projections (viewing angles) means 
that the PNS in the outer halo region is the most likely to 
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Figure 11. The 2D velocity fields for the face-on view (left) and 
for the edge-one one (right) for the model M3 at T = 4.5 Gyr. 
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same as Figure 11 but for the model M4. 



be observed as being strongly supported by rotation. Figure 
10 shows that a PNS with higher Kot/o" at i? = 2i?e is more 
likely to have higher Kot/o" at i? = 5iic: There is a corre- 
lation in stellar kinematics between the main body and the 
stellar halo of the elliptical, though the correlation is weak. 
Figure 10 also shows a stronger correlation between Kot /o"o 
at R = 2Rc and Vrot/o"o at R — 5Rc- The derived correla- 
tions imply that if elliptical galaxies are formed from major 
galaxy merging, the kinematic of their main bodies derived 
from integrated absorption spectra (for R < 2Rc) can be 
correlated with that of their outer stellar halos derived from 
radial velocity fields of PNe (for R > 5Rc)- The predicted 
"halo-host kinematic correlation" can be tested by future 
extensive systematic studies of PNe kinematics for elliptical 
galaxies with already known kinematics of the main bodies. 



3.2 Parameter dependences 

3.2.1 Orbital configurations 

The dependences of dynamical properties of PNSs in 
elliptical galaxy formed by major galaxy merging with m2 
= 1.0 on orbital configurations of the merging are described 
as follows. 
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Figure 13. The same as Figure 9 but for the different 5 merger 
models (M2, M3, M4, M5, and M6) viewed from three different 
directions (i.e., projected onto the X-Y, X-Z, and Y-Z plane). 



(1) Irrespective of orbital configurations, the radial 
profiles of the projected PNe number densities in PNSs can 
be fitted to the power-law profile for the entire halo regions 
of their host elliptical galaxies. The simulated PNSs of 
ellipticals in all models in the present study have the mean 
PNe densities more than an order of magnitude higher than 
those of stellar halos of their merger progenitor spirals. The 
central values of the projected PNe number densities do 
not depend on orbital configurations. These results suggest 
that structural properties of outer PNSs (thus stellar halos) 
can be significantly different between spirals and ellipticals. 

(2) Most of PNSs can show a significant amount of 
rotation in the outer halo regions [R > 5i?e) of their 
host elliptical galaxies, though the maximum values of the 
rotational velocities (Kot) depend on the viewing angles of 
the ellipticals. The 2D velocity fields of most PNSs shows 
minor-axis rotation in their halo regions, as seen in the 
main bodies {R < 2Re) of elliptical galaxies (e.g., Bendo 
& Barnes 2000). The details of the 2D velocity field of 
PNSs depend strongly on orbital configuration of galaxy 
merging. For example, the PNS of the elliptical formed 
from a retrograde-retrograde merging (model M4) shows 
a more inhomogeneous and more complicated 2D velocity 
field in its edge-on view than that from a prograde-prograde 
merging, though there are no significant differences in 
the 2D fields between the two for the face-on view (See 
Figures 11 and 12 and compares the two velocity fields 
with each other). The 2D velocity dispersion fields does not 
so strongly depend on orbital configurations than the 2D 
velocity ones. 

(3) Kot/o" of PNSs are more likely to be higher in 
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Figure 14. The same as Figure 10 but for the different 5 merger 
models (M2, M3, M4, M5, and M6) viewed from three diflferent 
directions (i.e., projected onto the X-Y, X-Z, and Y-Z planes). 



the outer halo regions of elliptical galaxies {R > 25 kpc 
corresponding to ~ 5i?o) than in the central regions _R ~ 5 
kpc, as shown in Figure 13. This result suggests that PNSs 
of elliptical galaxies axe likely to show rotational kinematics 
in their outer halo regions, if they are formed from major 
galaxy merging. Kot/o"o (rotation curve normalized to the 
central velocity dispersion cro) has the radial dependence 
similar to Kot/c (Figure 13). The flat rotation curve 
of Kot/fo seen in most of the models for their outer 
halo regions {R > 25 kpc ) suggests that such flattened 
shapes of Vrot/o"() is one of generic trends of stellar halos 
kinematics in elliptical galaxies formed from major merging. 

(4) There is a weak correlation between Kot /cr at R = 




R (kpc) 

Figure 15. The same as Figure 4 but for the three merger models 
with different (M2, M7, and M8). 



2Rc and Kot/o" at R — 5Rc in the sense that PNSs with 
larger Kot/o" at ii = 2Re are more likely to show larger 
Vroxjo at R = 5Re (See Figure 14). This result suggests 
that kinematics of outer stellar halos in elliptical galaxies 
can correlate with that of their main bodies, if elliptical 
galaxies are formed from major merging. Similar trend can 
be seen for Kot/co, which implies that elliptical galaxies 
with the main bodies rotating more rapidly are likely to 
show a larger amount of rotation in their outer halo regions. 
The derived two correlations, though weak, can be used as 
theoretical predictions that should be compared with future 
observations of PNe kinematics in elliptical galaxies. 

(5) These radial differences of Kot /cr characteristics of 
PNSs of elliptical galaxies formed by major merging (See 
Figures 13 and 14) can be derived only if we can obtain 
kinematical data for outer halo components (R ^ 5Re) of 
elliptical galaxies. This implies that we need to investigate 
observationally and theoretically kinematical properties of 
stellar halos beyond 2Rc to give some constraints on ellipti- 
cal galaxy formation and thus that PNe kinematical studies 
are ideal for this purpose. 



3.2.2 Mass ratios (1712) 

The dependences of dynamical properties of PNSs in merger 
remnants on are described as follows. 

(1) The projected PNe number densities ppNo 
(logio(N/kpc^)) are globally lower in the remnants with 
smaller m2, in particular, for the outer parts {R > 5Re 
corresponding to ~ 25 kpc) of PNSs (See Figure 15). This 
means that more flattened elliptical galaxies are more likely 
to show lower ppNc, because the morphological properties of 
the host galaxies formed from these mergers with lower rn2 
show more flattened intrinsic shapes. The power-law slope 
of the density proflle ppNe (logio(N/kpc^) does not strongly 
depend on 7712 for R < 5Re, however, there is a clear sign of 
steeper slope in the model with m2 — 0.1 for R > 5-Rc owing 
to the rather low PNe density there. The result implies that 
the ppNe is steeper in Es formed from major galaxy merg- 
ing (e.g., m2 > 0.3) than in SOs from unequal-mass/minor 
galaxy merging with (0.1 < m2 < 0.3). 

(2) The kinematics of PNSs in E/SOs formed from 
galaxy merging with lower m2 show more rapid rotation, 
larger maximum values of Kot , and steeper radial gradients 
of Kot for R < 2Re (See Figure 16). Given the fact that 
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Figure 16. The radial profiles of Vl-ot/o" (upper) and the normal- 
ized rotation curve Vrot/cro (lower) for the unequal-mass merger 
model M8 with m2 = 0.3 viewed from three different viewing 
angles (projected onto the X-Y, X-Z, and Y-Z planes). 



remnants of galaxy merging with smaller 7712 show more flat- 
tened shapes, this result suggests that more flattened E/SO 
have higher Kot and steeper radial Kot gradients. Irrespec- 
tive of m2, the outer regions of PNSs {R ^ 5-Rc) can show 
high V^ot and thus large Vrot/o", owing to the nearly flat ro- 
tation curve and the radially decreasing a (See Figure 17). 
This suggests that stellar halos in E/SOs formed by galaxy 
merging show different kinematics compared with that of 
the Galaxy with little rotation. 

(3) The larger values of Vrot /cr in the PNS for edge-on 
projections (i.e., X-Z and Y-Z) in SOs formed from unequal- 
mass merging suggests that rotational terms in Jeans equa- 
tions should be considered when we estimate total masses 
of SOs using PNe data (See Figure 16). The small values of 
Kot/o" for the face-on projection (i.e., X-Y) imply that the 
total mass of an SO (thus mass-to-light-ratio) can be signif- 
icantly underestimated if it is viewed from face-on and if 
the mass is estimated by the Jeans equation with correc- 
tion terms of rotation (and if the projected central velocity 
dispersion does not depend on the viewing angle). The ob- 
served small M/L in some early- type galaxies (e.g., in NGC 
3379) could be due to this viewing angle effect. 



4 COMPARISON WITH OBSERVATIONS OF 
NGC 5128 PNE SYSTEMS 

4.1 Selection of the best model 

We present the results of one of the best models (model 
M9) in which the observed structural and kinematical prop- 
erties of the PNS in NGC 5128 (PFF04a) are the most self- 
consistently reproduced. We select this best model among 
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Figure 17. The same as Figure 1 but for the model M9. opNc 
of 28.2 PNe per Lq in S-band is adopted so that the PNe num- 
ber density at the effective radius of the merger remnant can be 
consistent with the observed one for the NGC 5128 PNS. 
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Figure 18. The same as Figure 3 but for and the observational 
data (thin) by PFF04a and the model M9 (thick). 



the 32 models investigated in the present study by checking 
whether the following five fundamental observational results 
(e.g., Israel 1998; PFF04a) can be reproduced reasonably 
well in each model: (1) Rs = 5.2 kpc. Mi (total luminous 
mass) = 1.4 x lO'^^ Mq (for Mi/Lb = 3.5), and spherical 
appearance of the elliptical galaxy, (2) the projected PNe 
number density ppNe (logio(N/kpc'^) ~ 1.0 at 2Rc and the 
slope of the power-law profile similar to —2.5, (3) the rota- 
tion curve that rises till R ~ 2Re and becomes flat at ii > 
2Re with the maximum Kot of ~ 100 ± 20 km s~^, (4) the 
central velocity dispersion ctq of ~ 140 ± 10 km s"'^, and (5) 
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Figure 19. The same as Figure 5 but for the observational data 
(left) and the model M9 (right). 



the 2D velocity field showing both minor and major axes ro- 
tation with the zero velocity curve/contour (ZVC) twisting 
significantly. 

Owing to the purely coUisionless nature of the present 
simulations, we can match any model to the above observa- 
tion (1) by rescaling the size and the mass of the simulation 
Therefore, the above (2) - (5) constraints can be used in se- 
lecting the best possible model in the present study. Major 
merger models with m2 = 1.0 can explain the above (2), 
(4), and (5) reasonable well, however they have difficulties 
in explaining (3) owing to slowly rising Kot with smaller 
maximum values of V^ot ■ Unequal-mass merger models with 
m2 = 0.1 — 0.3 can become early-type E/SOs, however, they 
are very fiattened in shapes, strongly supported by rotation, 
and have 2D velocity fields that are not similar to the ob- 
served ones. These models accordingly can not explain the 
above (1), (4), and (5) in a fully self- consistent manner axiA 
thus can not be regarded as the best model. Thus one of the 
best models can be unequal-mass mergers with m2 of some- 
where between 0.3 and 1.0 (viewed from a certain direction): 
However, it should be stressed here that we possibly could 
miss out the major merger model with m2 = 1.0 that can 
explain the above five points self-consistently, owing to the 
limited number of the simulated models. 

The best model for which we show the results below 
is the model 9 with m2 = 0.5 in which a highly inclined 
orbital configuration and a smaller pericentre distance are 
assumed. We investigate dynamical properties of the PNS in 
this model at T = 4.5 Gyr when the remnant is dynamically 
well relaxed to show regular distributions of PNe both in the 
2D distribution of the projected number density and in the 
radial one (See Figures 17 and 18). This model with smaller 
pericenter distance and a highly inclined disk with respect 
to the orbital plane can form stellar shells and gaseous rings 
perpendicular to the major axis of the merger remnant if 
gaseous dissipation is included (Bekki 1997, 1998a). Thus, 
our future more sophisticated model with gaseous dissipa- 
tion and star formation will be capable of explaining the 
observed HI distribution perpendicular to the photometric 
major axis and the outer shells as well as the above five ob- 
servations. Dynamical properties of globular cluster systems 
and fine structures observed in NGC 5128 (Peng et al. 2002; 
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Figure 20. A comparison between the observed ZVC (zero veloc- 
ity curve) and the simulated one for the model M9. The observed 
one (PFF04a) is represented by a solid line and the filled circles 
indicates the location of cells with |Vios| (the absolute magnitude 
of lino-of-sight- velocities) ^ 3.0 km . The locations of the cells 
with X < kpc and Z < kpc and those with X > kpc and Z 
> kpc are not plotted for the simulation so that we can more 
readily infer the possible ZVC of the simulation model from the 
distribution of the cells with |Vios| ^ 3.0 km . The mathemat- 
ical expression of the observed ZVC is given in the main text. 
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Figure 21. The same as Figure 6 but for the observational data 
(left) and the model M9 (right). 



2004b, c) will be discussed in the best model(s) of our future 
studies (Bekki & Peng 2005, in preparation). 



4.2 Advantages and disadvantages of the best 
model 

Figures 19 presents the comparison between the simulated 
2D velocity field of the PNS of the best model and the ob- 
served one. In order to match the observation with the sim- 
ulation, the observed coordinate —X {Y) for the PNS in 
NGC 5128 (PFF04a) is set to represent +X [Z) in Figure 
19 (and 20, 21) for convenience. The best model is consis- 
tent with the observation at least qualitatively in the sense 
that it shows (1) strong major-axis rotation extending to 35 
kpc (~ 7-Ro), (2) weaker but significant rotation along the 
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Figure 22. The same as Figure 7 but for the observational data 
(thin) and for the model M9 (thick). In order to compare the 
simulation result with the observations in a more self-consistent 
manner, the way of binning is exactly the same between the sim- 
ulation and the observation. 



minor axis, and (3) the line of zero velocity (referred to as 
zero velocity curve, ZVC) both misaligned and twisted with 
respect to the major axis of the PNe distribution. 

Figure 20 more clearly describes how the ZVC looks 
like in the 2D velocity field and whether it is consistent 
quantitatively with the observed ZVC. The observed ZVC 
can be parameterized as follows (PFF04a): 

-J!:(0.125|X| -F5) 



Z = 



+ 2.56) 



(2) 



where X and Z arc exactly the same as the coordinate X and 
Z in the simulations. As shown in Figure 20, the simulated 
ZVC starts twisting at Z ~ ± 10 kpc whereas the observed 
ZVC starts twisting at Z ~ ± 4 kpc. The direction of the 
simulated ZVC for \X\ > 10 kpc is broadly consistent with 
the observed one, though the simulated ZVC is more aligned 
with the minor axis of the 2D PNe distribution (shown in 
Figure 19) compared with the observed one for |X| ^ 10 
kpc. Those less successful reproduction of the model sug- 
gests that (1) some physics which are not included in the 
present study (e.g., gaseous dissipation and star formation) 
should be considered to reproduce fully self-consistently the 
observations and (2) we need to explore wider sets of model 
parameters (e.g., merger orbits) to find a fully self-consistent 
model. 

Figure 21 presents the comparison between the simu- 
lated 2D velocity dispersion field of the PNS and the ob- 
served one (PFF04a). The simulation is consistent broadly 
with the observation in the sense that (1) it shows an inner 
flattened shape in the 2D velocity dispersion field and (2) 
radial gradient of velocity dispersion along the minor axis 
is steeper than that along the major axis. Interestingly, the 



observation shows an isolated region with high velocity dis- 
persion (~ 200 km s~^) around {X, Z) = (15, -5) (kpc). 
An apparently isolated region with moderately high velocity 
dispersion (~ 120 km s~^) can be discernibly seen also in 
the simulation around (X, Z) = (25, —2) (kpc). Although 
it is not so clear whether the presence of such local, dynam- 
ically hot regions in the 2D velocity dispersion fields have 
some physical meanings of elliptical galaxy formation, more 
extensive comparison between observations and simulations 
in terms of the locations of the isolated, dynamically hot 
regions should be made to determine the best model after 
a larger set of PNe kinematical data become available, in 
particular, for the regions with X > kpc and Z < —IQ 
kpc. 

Figure 22 demonstrates that the best model can ex- 
plain the observed rotation curve (Kot) profile reasonably 
well for J? < 3 kpc and 10 < J?50 kpc. However there is 
a significant difference in Kot profiles between the observa- 
tion and the simulation for R ^ 5 kpc: The observational 
rotation curve rises more rapidly than the simulated one for 
the inner region of NGC 5128. The radial profile of veloc- 
ity dispersion (u) in the model is systematically higher (a 
factor of 1.2 at _R ~ 2Rc) than the observed one. Wo sug- 
gest that these less successful reproduction of the best model 
for Kot around R = 5 kpc and the a profile is due to the 
model's not including PNe formation from gas: All PNe are 
assumed to originate from coUisionless old stellar disks. We 
can expect that if PNe can form from gas, new PNe have a 
larger amount of rotation and smaller amount of random ki- 
netic energy owing to efficient energy dissipation of random 
kinematic energy during galaxy merging. We thus suggest 
that our future more sophisticated model with gaseous dis- 
sipation can more successfully reproduce the observed radial 
profiles of Kot and a. 



5 DISCUSSION 

5.1 Constraints on galaxy formation from stellar 
halo kinematics 

PFF04a have revealed that the PNS of NGC 5128 has a sig- 
nificant amount of rotation with V/a between 1 to 1.5 in its 
outer stellar halo region {R > 5Re). The present study has 
shown that outer stellar halos in most of major merger mod- 
els have a significant amount of rotation and thus suggested 
that a rotating stellar halo seen in NGC 5128 is not just a 
special case but a rather general trend of elliptical galaxies, 
if ellipticals are formed by galaxy merging. The kinematics 
of PNSs in several elliptical galaxies have been investigated 
so far (CiarduUo ot al. 1993; Hui et al. 1995; Arnaboldi et al. 
1998; Mendez et al. 2001; Romanowsky et al. 2003; Napoli- 
tano et al. 2004; PFF04a) and found to show rotation in 
some of the ellipticals. Although these observational studies 
help the authors to provide a reasonable dynamical model 
(e.g., total mass of a galaxy) for each individual case, the 
total number of PNSs investigated is too small for them to 
make any robust conclusions on the general trend of kine- 
matics in outer stellar halos (R > 5Re) of elliptical galaxies. 

We suggest that extensive systematic studies of kine- 
matical studies of PNe for 2Re ^ 71 ^ 10-Re for a larger 
number of elliptical galaxies are doubtlessly worthwhile, be- 
cause the kinematics of outer stellar halos can provide strong 
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Table 2. Mass comparison 



Model no. Fitted mass (Mt.f) Actual mass (Mt, a) A^t,f/A^t,A 

Ml (edge-on) 1.7 x IO^Mq 8.6x10"Mq' 0.20 

Ml (face-on) 7.4 x IO^OMq 8.6 x 10" Mq 0.09 

M8 (edge-on) 9.3 x 10"Mq 7.1 x 1O"M0 1.31 

MS (face-on) 1.5 x lO^^Mo 7.1 x 10"M© 0.02 



constraints on elliptical galaxy formation models, such as 
the monolithic collapse scenario (e.g., Larson 1974; Carlberg 
1984) the merger one (e.g., Toomre 1977; Barnes 1992), and 
multiple merging of subgalactic clumps or dwarfs (e.g.. Cote 
et al. 2002). Although kinematical properties of outer stellar 
halo have not been investigated by the monolithic collapse 
scenario (Larson 1974), it is unlikely that the outer stel- 
lar halos in elliptical galaxies formed in monolithic collapse 
show more rapid rotation than the inner main bodies, be- 
cause random kinematical energy is likely to be dissipated 
away more efficiently to form a stellar system with rotation 
in the inner regions of protogalaxies where gas densities are 
higher (thus more dissipation is highly likely). If elliptical 
galaxies are built up through multiple merging of dwarfs 
with random orientation of merging, there appears to be 
no strong physical reasons for them to show a significant 
amount of rotation in their outer stellar halos. Thus stellar 
halo kinematics that can be proved by kinematical studies 
of PNSs in currently ongoing observations can be used to 
assess the viability of each of the above three scenarios of 
elliptical galaxy formation. 

5.2 Halo-host connection 

The present study discovered that both the radial profiles of 
projected PNe number densities (ppNe) and the kinematical 
properties of PNSs in merger remnants depend strongly on 
the mass ratios of the merger progerntor spirals (See Figures 
14, 15, and 16). Given the fact that morphological properties 
of merger remnants are demonstrated by numerical studies 
to depend strongly on the mass ratios (e.g., Bekki 1998b; 
Naab & Burkert 2003) , the above discovery can provide the 
following three predictions on the possible correlations in 
dynamical properties between PNSs and their host galax;- 
ies. Firstly, more flattened E/SO galaxies are more likely to 
have smaller mean ppNc and higher K-ot/o" at 5Rc, because 
unequal-mass galaxy mergers with rn2 — 0.3 — 0.5 finally 
become more flattened E/SOs (Bekki 1998b; Naab & Burk- 
ert 2003) . Secondly, stellar halos of SOs with thick disks are 
more likely to show a significantly smaller ppNe compared 
with those of Es, because SOs with thick disks can be formed 
from lower mass-ratio unequal-mass galaxy mergers (m2 — 
0.1 — 0.3). Thirdly, E/SO galaxies can show a correlation in 
kinematics between V/a (e.g., at ~ 2Re) of the main bod- 
ies and that of their stellar halos (e.g., at ~ 5Re), if most 
of them arc formed from galaxy merging. 

These three predictions will be more readily tested in 
future and ongoing kinematical studies of main bodies and 
stellar halos of E/SO galaxies by the planetary nebula spec- 
trograph (Douglas et al. 2002) and integrated field spec- 
trographs like SAURON (Bacon et al. 2001). Furthermore, 
as have been demonstrated by HHOO an PFF04a, both the 
MDF and the kinematics of the stellar halo in NGC 5128 
can be significantly different from those of the stellar halo 



of the Galaxy in the sense that the stellar halo of NGC 

5128 are more metal-rich and more strongly rotating. These 
differences in stellar halo properties may well be remark- 
able and fundamental differences between spiral and ellipti- 
cal gala^xies, though extensive observational studies of PNSs 
in spirals beyond the Local Group that can prove dynamical 
properties of the stellar halos have not been yet conducted. 
Although the MDF was investigated for the stellar halo of 
the edge-on SO NGC 3115 (Elson 1997; Kundu & Whitmore 
1998), systematic studies on the MDFs of stellar halos and 
dynamical properties of PNSs in SOs have not been done yet. 
Thus future observational studies of structural and kinemat- 
ical properties of PNSs in spirals and SOs, combined with 
those ongoing for ellipticals, will provide a new clue to the 
origin of the Hubble sequence. 

5.3 Mass Estimates 

One ongoing field of study using PNSs is the total mass of 
elliptical galaxies. Typically, analysis of this sort use the in- 
tegrated rotation and velocity dispersion profiles of the PNe 
to derive a total meiss within a given radius (e.g. PFF04a, 
Romanowsky et al. 2003, Mondez et al. 2001). As one can see 
from Figure 9, a single galaxy can have a wide range of rota- 
tion and dispersion profiles depending on the viewing angle, 
and this could plausibly lead to biases in the mass estima- 
tion. We test this by adopting the mass estimation procedure 
outlined in PFF04a, which uses the spherical Jeans equation 
and assumes an isotropic distribution of orbits. While this 
is a simple assumption, it is one that is typically used and 
there has yet to be any strong kinematic evidence that or- 
bits in ellipticals are strongly anisotropic. When we fit the 
kinematic profiles of the major merger model from the face- 
on and edge-on views, we find that the face-on view gives 
a mass that is a factor of two lower than the edge-on view. 
The reason for this can be seen qualitatively in Figures 7 
and 8 which shows that the rotation and velocity dispersion 
profiles are both lower for the face-one view. 

Even more intriguing, the simulations consistently show 
that the dark matter is at a higher velocity dispersion than 
the stars. In fact, while the face-on view under-represents 
the mass relative to the edge-on view, both viewing angles 
underestimate the total mass. In order to estimate the total 
mass of a merger remnant from PNe kinematics and thereby 
discuss this point in a more quantitatively, we use the same 
method as those used by Hui et al. (1995) and PFF04a: We 
apply the spherical Jeans equation to the major-axis rota- 
tion and line-of-sight velocity dispersion profile (e.g.. Figures 
7 and 8) in a merger remnant (i.e., an elliptical galaxy) to 
derive the dynamical mass of the remnant (See PFF04a for 
more details on the methods of mass estimation from PNe 
kinematics). 

Table 2 summarizes the actual masses of the simu- 
lated E/SOs within 45 kpc (Mt,a) and the masses estimated 
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from the simulated radial velocity dispersion profiles and 
the Jeans equation used in Peng et al. (2004a) within 45 
kpc (Mt,f) for two interesting cases. Compared to the total 
amount of mass in the simulation, our simple mass estima- 
tion from the stellar data viewed edge-on (best case) un- 
derestimates the total mass within 45 kpc by a factor of 5. 
In one model for a minor/unequal-mass merger viewed face- 
on, the total mass is underestimated by a factor of 50 and 
requires no dark halo. If ellipticals result from these types 
of mergers, then this could be a very important considera- 
tion, especially considering that current mass estimates for 
ellipticals axe surprisingly low. This matter warrants further 
investigation with a variety of different elliptical formation 
models. (Dckcl ct al. 2005 have recently discussed this prob- 
lem with a complementary approach and also find masses 
can be underestimated). However, we point out that the 
supposedly low M/L recently measured are for intermediate- 
luminosity galaxies (e.g. NGC 3379 and NGC 5128) and not 
for all galaxies. If this "stellar kinematic bias" is the reason 
for anomalously low mass estimates, then it may be more 
pronounced in lower luminosity ellipticals, providing a clue 
to their formation histories. 



6 CONCLUSION 

We have numerically investigated structural and kinemati- 
cal properties of PNSs in elliptical galaxies formed by galaxy 
merging in order to elucidate the origin of the observed 
physical properties of PNSs. We have mainly investigated 
the radial profiles of projected number densities, rotational 
velocities, and velocity dispersion of PNSs and the two- 
dimensional velocity fields of PNSs in elliptical galajdes. We 
also compared the simulated dynamical properties of PNSs 
with the corresponding observations of the PNS in NGC 
5128 and thereby tried to provide the best merger model for 
the PNS in NGC 5128. We summarize our principal results 
as follows. 

(1) The radial densities profiles of PNSs can be fitted to 
the power-law in the entire halo regions of elliptical galaxies 
formed from major merging (m2 = 1.0). The projected PNe 
number densities (ppNe) at -R ~ 5Re in elliptical galaxies 
are more than an order of magnitude higher than those of 
the halos of the merger progenitor spirals. These results 
do not depend strongly on the model parameters of major 
merging. The main reason for rather high ppNe of PNSs is 
that a significant fraction 10%) of disk stars are stripped 
from the disks and redistributed in the halo regions during 
major merging. These results suggest that elliptical galaxies 
have higher PNe number densities (ppNe) in their halos 
than spiral ones. 

(2) PNSs of elliptical galaxies formed from major 
merging can show a significant amount of rotation (Kot/o" 
> 0.5) in their outer halo regions {R > 5Rc)- The derived 
rotation in the outer halos results from the angular mo- 
mentum redistribution of disk stars during galaxy merging 
(i.e., conversion of orbital angular momentum of merging 
two spirals into the intrinsic one of the merger remnant). 
Kot/o" of PNSs in ellipticals are more likely to be larger 
in the outer parts than in the inner ones, though radial 
profiles of Kot/cr are diverse between different models. 



Vrot /cr at 5Rc can weakly correlate with Kot/o" at 2Rc 
for the PNSs in such a way that PNSs with larger Kot/o" 
at 2Re are likely to show larger Kot/c at 5Re. This 
result implies that the kinematics of outer stellar halos in 
elliptical galaxies can correlate with that of the main bodies. 



(3) The two-dimensional (2D) fields of velocity and 
velocity dispersion of PNSs in elliptical galEixies formed by 
major merging are quite diverse depending on the orbital 
configurations of galaxy merging, the mass ratios of the 
progenitor spirals, and the viewing angles of the galaxies. 
For example, PNSs of ellipticals formed from prograde- 
prograde merging can show more rapid rotation along the 
major axes of the PNe density profiles in their 2D velocity 
fields compared with the retrograde-retrograde ones. For 
most of the models, the 2D velocity fields show minor axis 
rotation in the outer halo regions of elliptical galaxies. It is 
doubtlessly worthwhile for future observational studies to 
systematically investigate dynamical properties of PNS for 
a larger number of elliptical galaxies and thereby to confirm 
whether the predicted uniformity and diversity in the 2D 
velocity fields of PNSs can be seen in elliptical galaxies. 



(4) ppN of PNSs are more likely to be lower in E/SOs 
formed from galaxy merging with smaller mass ratios (e.g., 
1712 = 0.3, i.e., unequal-mass merging). Furthermore, PNSs 
in the remnants of mergers with smaller m2 have more 
flattened shapes and a larger amount of rotation (Kot) both 
in their main bodies and in their outer halos. Therefore 
these results suggest that more flattened E/SOs are more 
likely to show both lower ppNe and higher Vrot (for a 
given luminosity range of E/SOs). This predicted correla- 
tion can be readily confirmed in future observations of PNSs. 



(5) The observed kinematical properties of the PNS 
in NGC 5128 (e.g., Vrot/cr between 1 and 1.5 and minor 
axis rotation) can be broadly consistent with the present 
best model with m2 = 0.5, a small impact parameter, and 
highly inclined initial disks (similar to a polar-orbit). The 
observed "kink" in the zero velocity curve (ZVC) of the 2D 
velocity field of the PNS in NGC 5128 can be reproduced 
reasonably well, though the location where ZVC begins to 
twist is different between the best simulation model and the 
observation. Some disadvantages of the present best model 
in explaining self-consistently the observed kinematics of 
the PNS suggest that gas dynamics and star formation may 
well play an important role in the formation of the PNS 
observed in NGC 5128. 



(6) The mass estimates of merger remnants viewed 
face-on are likely to be a factor of two lower than those 
viewed edge-on. However, even the mass estimates of sys- 
tems viewed edge-on can be low by a factor of 5, and in the 
worst case (face-on E/SO model) can be low by a factor of 
50. If this conclusion is applicable to real early- type galax- 
ies, it has interesting consequences for current observational 
work. 
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APPENDIX A: DERIVATION OF 2D DENSITY 
AND VELOCITY FIELDS 

In order to compare the simulated 2D density and velocity 
fields with the observed ones in a more self-consistent man- 
ner, we adopt the methods that are quite similar to those 
used for kinematical analysis of observational data of PNSs 
in elliptical galajcies (e.g., PFF04a). We first determine the 
position angle (^p) of the major axis of the stellar mass dis- 
tribution projected onto the x-z plane (i.e., perpendicular to 
orbital plane of galaxy merging) for each model at final time 
step. Then we rotate the model by —6p within the x-z plane 
so that the major axis of the mass distribution can coincide 
with the X-axis. For convenience, the coordinate {x, y, z) is 
renamed as [X, Y , Z) after the rotation of the model. We 
represent the kinematics of a model by estimating line-of- 
sight (LOS) velocity moments as a function of position with 
a nonparametric smoothing algorithm. At the position of 
each stellar particle, we apply a local linear smoother using 
a Gaussian kernel function with the smoothing length of 0.17 
in our units (corresponding to 3 kpc). We choose this value 
of 3 kpc, firstly because recent observational study (PFF04a) 
chose this value and secondly because we can more clearly 
see the detail of the 2D velocity field of a model throughout 
the outer halo region (extending to i? ~ 10J?e without losing 
resolution within Re. 

We divide the entire halo region with the size of ARa 
(= 70 kpc corresponding to ^ 14_Rc) into 50 x 50 cells for a 
merger model in each projection and estimate linc-of-sight- 
velocity of v{Xi,Zi) for the X — Z projection {v{Xi,Yi) 
and v{Yi,Zi), for {Xi,Yi) and (Yi,Zi), respectively), at the 
center of each cell, i.e., {Xi,Zi) (i = 1 — 50) based on the 
smoothed velocity filed (described above). Total number of 
cells in each projection is fixed at 2500 for all merger mod- 
els in the present study, because we can more clearly see the 
global changes of 2D velocity field without suffering small- 
scale, and rapid variation resulting from small number of 
particles in each cell for this cell number. We mainly show 
the 2D velocity fields projected onto the X-Y plane (referred 
to as the face-on view for convenience) and the X-Z plane 
(the edge-on plane). The same method is used for estimation 
the 2D density field of the PNS in model. 

In converting the 2D stellar mass density field into the 
2D PNe one, we assume that the luminosity-specific PNe 
(number) density represented as apNe is 9.4x10"^ PNe/L© 
(in B-band) This represents the rmmber of PNe expected 
in the brightest 2.5 magnitudes of the PNLF, which is also 
roughly equivalent to the typical PN survey depth. Observa- 
tions of MSl's PNS (Ciardullo et al. 1989) showed that apNe 
can range from 2.9xl0~^to 39.3x10"^, and accordingly the 
adopted value above is consistent with these observations. 
This value of 9.4xl0~®, which originates from the 2nd col- 
umn of the table 4 in Ciardullo et al (1989), is used for apNe 
throughout this paper unless specified. While this is very 
obviously a simplification — a is known to vary as function 
of metallicity, and observational incompleteness will often 
vary across a galaxy — it is a good first-order approximation 
to the true PN population for our purposes. We also assume 
that M/Lb = 3.5 for stellar components of galaxies in all 
models. 

The total number of PNe (iVpNe) within 5Re of the sim- 
ulated elliptical galaxies is typically ~ 350. As shown in Fig- 



ure 1, even for the small number of star particles (1000) in 
the halo of the disk, the initial spherical distribution of PNe 
(stellar halo) can be reproduced reasonably well. In most of 
the merger models, typically 10 % of the total stellar mass 
in the merger remnant can distribute throughout the halo 
regions {R > 2Re). Therefore, more than 10* stellar particles 
(up to 5 X 10'') can be used to derive the smoothed density 
(velocity) field for a model. Thus we can derive the smoothly 
changing density (velocity) fields throughout the entire halo 
regions from original, more discrete density (velocity) ones 
for the adopted total number of stellar particles in simula- 
tions by using the smoothing method described above. 

We also derive the radial profiles of rotational velocity 
(Kot) and velocity dispersion (a) for each model based on 
the same method as used for observational analysis on the 
NGC 5128 PNS (PFF04a). For example, when we try to de- 
rive the radial profile of Kot along the major axis of a PNS , 
we use PNe within either a perpendicular distance of ± 2 kpc 
from the major axis or a ±10° cone centered on the major 
axis. This method is adopted so that more PNe can be in- 
cluded in the analysis of the radial profiles for the outer halo 
regions, where less number of stellar particles can prevent us 
from making a reasonable estimate of Kot and a. The errors 
in Kot (o") shown in figures of the present paper (e.g.. Figure 
7) are equal to Kot/-\/2(iV - 1) (<T/^2(iV - 1), where N is 
the total number of particles for a given radial bin. We do not 
intend to separately estimate Kot and a for PNe originating 
from spiral's stellar halos (thus metal-poor stars) and those 
from disks (metal-rich ones), partly because observational 
data sets on metallicities of PNe are not currently available. 
Although we find that there are no significant kinematical 
differences between metal-poor and metal-rich PNe in some 
models, we will discuss possible metallicity dependences of 
PNe kinematics in our forthcoming papers. 

In order to investigate the smoothed 2D velocity fields, 
we use the same smoothing method as that used for deriv- 
ing the 2D density fields of PNSs. The method of velocity 
smoothing is described as follows. The particle with the lo- 
cation X, and the velocity V, and the mass of M is con- 
sidered to be composed of Ns particles (hereafter referred 
to as "smoothing particles") that are located within R (cor- 
responding to smoothing length) from the particle ( "parent 
particle") and have position vectors of x, (i = 1 ~ iVs) 
with respect to the parent, velocity vectors of V (i.e., the 
same as that of the parent), and masses of M/Ns. The spa- 
tial distribution of these smoothing particles with respect to 
the parent particle (X) follows the Gaussian distributions. 
Therefore the probability of a i-th smoothing particle with 
the position of X + Xi is proportional to exp( — ri'^/2i?^), 
where is the distance between the parent particle and the 
smoothing one. If the number of smoothing particles for a 
stellar particle is 100, the total number of "particles" used 
for 2500 bins is about 10®. Therefore, each bin includes 400 
"particles" on average. 

For each bin, we average the velocities of "particles" 
with the positions within the bin and derive the velocity at 
the location of the bin. The adopted smoothing length is con- 
sistent with that used for observations by PFF04a. Also we 
confirm that this method can reproduce the "spider shape" 
(which is characteristic of rotational kinematics) of the 2D 
velocity field of the initial disk. Thanks to this smoothing 
method, we can derive 2D density, velocity, and dispersion 
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fields of PNSs that can be tested against tlie corresponding 
observations in a fully self-consistent manner. 

This study is based totally on dissipationless simula- 
tions of galaxy mergers so that it can not discuss struc- 
ture and kinematics of PNc formed from gas during/after 
galax;y merging. Accordingly, the present results can be more 
reasonably compared with the observed E/SOs formed from 
"dry" (dissipationless) mergers that do not create new PNe. 
In our future papers, we will discuss how the introduction 
of PNe formation in our numerical simulations can change 
main conclusions derived in the present study. 



els are (1) shallower slopes of radial a profiles in multiple 
merger models and (2) larger a differences between dark 
matter halo and PNe for the multiple merger model. The 
above (1) suggests that there can be diversity in radial a 
profiles in elliptical galaxies, if most of elliptical galaxies are 
formed either from pair mergers or from multiple ones. The 
above (2) implies that radial a profiles of PNSs alone can 
not allow us to correctly estimate of the total masses of 
elliptical galaxies formed by multiple merging. We have a 
plan to investigate this problem related to mass estimation 
of galaxies by PNe kinematics in a more extensive manner 
in our forthcoming papers (Bekki & Peng 2005) 



APPENDIX B: MULTIPLE MERGERS 

Each multiple merger model contains equal-mass spiral 
galaxies with random orientations of intrinsic spin vectors 
which are uniformly distributed within a sphere of size 
6-Rd. The most important parameter in this multiple merger 
model is the ratio of the initial kinematic energy (Tkin) of 
the merger to that of initial potential (W). By varying this 
ratio (tv; defined as |2rkin/H^|) from 0.25 to 0.75, we in- 
vestigate how tv controls the final PNc kinematics of the 
merger remnants. We mainly investigate the results of two 
extreme cases: (1) where the initial kinetic energy of a mul- 
tiple merger is due entirely to the random motion of the five 
constituent galaxies (referred to as "dispersion supported") 
and (2) where it is due entirely to (rigid) rotational motion 
("rotation supported"). Our investigation of these two cases 
enables us to understand how the initial rotation (disper- 
sion) can control the final kinematical properties of PNe in 
merger remnants. 

Although wc have derived the results of three models 
with tv = 0.25, 0.5, and 0.75 for each case (i.e., six mod- 
els in total), we describe the result of the "rotation sup- 
ported" model with tv = 0.5 (labeled as MIO in the Tar- 
ble 1). This is firstly because this multiple merger model 
show some interesting differences in PNe kinematics com- 
pared with pair merger models, and secondary because this 
model shows typical behaviors in PNe kinematics among 
multiple merger models. The merger remnant of this rota- 
tion supported model shows a flattened stellar (thus PNe) 
distribution if it is seen from edge-on and have effective ra^ 
dius of 11.4 kpc for the stars and 29.8 kpc for the dark 
matter halo. 

We briefly summarize the results as follows. The radial 
gradient in a is significantly shallower in the multiple merger 
models than in pair merger models. These shallower a pro- 
files in PNe kinematics can be seen in most multiple models 
and thus can be regarded as one of characteristics of PNSs 
in multiple merger remnants. The difference in a at each ra- 
dial bin between PNe and dark matter halo is slightly larger 
than that seen in pair merger models Both dark matter halo 
and PNe show larger rotational velocities in the outer part 
of the halo and PNe can show larger Kot than dark matter 
halo. 

Thus PNSs in elliptical galaxies formed by multiple 
mergers show similar kinematics to those in elliptical galax- 
ies by pair mergers, though the inaxiinum values of a and 
Kot can be different between the two cases owing to the 
larger masses in multiple mergers. Remarkable differences 
in PNe kinematics between the two different merger mod- 



